Mathematical Biosciences and Engineering

Following the suggestion of L. F. Jaffe [1] we propose a mathematical model of fast calcium induced calcium influx waves (CICI Waves). They can propagate at relatively high speeds (up to 1300 micrometers/s). According to [1], they propagate due to a mechanochemical interaction of actomyosin network with the cell membrane. The local stretching of the membrane caused by actin filaments opens mechanically operated ion channels resulting in the influx of calcium to the cell. Moreover, stretching a cells membrane at one point opens nearby stretch activated calcium channels because the mechanical force is relayed by the actin filaments interconnected by myosin bridges. The number of bridges as well as filament density increases with calcium concentration, causing the contraction of the actomyosin network. Thus, the force acting on the membrane from tangled actin filaments is transmitted ahead of the moving front of the calcium concentration. As a result, the ion channels are opened even before the signal of calcium reaches them. This leads to much larger propagation speed of CICI waves in comparison with calcium induced calcium released (CICR) waves, where the wave is sustained by the diffusion of calcium and autocatalytic release of calcium from the internal stores (e.g. endoplasmic reticula).

Evolutionary cancer therapy (ECT) applies principles of evolutionary game theory to prolong the effectiveness of cancer treatment by curbing the development of treatment resistance. It was shown to increase time to progression while decreasing the cumulative drug dose. ECT individually tailors treatment schedules for patients based on their cancer dynamics and, thus, requires regular follow-up and precise measurements of the cancer burden. The current literature on ECT often overlooks clinical realities, such as rather long intervals between tests, possible appointment delays and measurement errors, in the development of the treatment protocols. In this study, we assess the clinical feasibility of ECT for metastatic non-small cell lung cancer (NSCLC). We create virtual patients with cancer dynamics described by the polymorphic Gompertzian model, based on data from the START-TKI clinical trial. We assess the effects of longer test intervals, measurement error and appointment delays on the expected time to progression under the evolutionary therapy protocols. We show that a higher containment level, although it increases time to progression in the models predictions, may lead to premature treatment failure in the presence of measurement error and appointment delay. Further, we show that the ECT protocol with a single containment bound is more robust to the clinical realities than the protocol with two bounds. Finally, we show that a dynamically adjusted treatment protocol can be beneficial for individual patients, but requires a thorough follow-up. This study contributes to the design of a clinical trial and the future clinical implementation of evolutionary therapy for NSCLC.

The left ventricle (LV) exhibits torsional deformation during systole, and mechanical relaxation begins during the isovolumic phase. Recent advances in imaging techniques, such as MRI, have revealed that myocardial tissue deformation and sarcomere length changes occur during the isovolumic relaxation phase, even when the chamber volume remains constant. Although such ventricular deformation during the isovolumic phase is considered important for blood ejection and filling efficiency, its mechanistic contribution to contraction and relaxation remains unresolved. In this study, we hypothesized that sarcomere length dynamics during the isovolumic phase affect the isovolumic contraction and relaxation time (IVCT and IVRT) by regulating the contraction force via the force-velocity relationship of ventricular myocytes. To investigate this hypothesis, we focused on experimentally reported differences in the relationship between sarcomere length and LV volume across the endocardial and epicardial layers, as described by Rodriguez et al. We constructed and compared two types of hemodynamic models within the same integrated framework consisting of a circulation model, a LV model, and a myocardial cell contraction model by Negroni-Lascano et al., which differ only in how sarcomere length is determined: a volume-based length model (VL model), in which sarcomere length is uniquely determined by LV volume, and a volume-force-coupled length model (VFL model), in which sarcomere length is determined by the balance between LV volume and contraction force. Simulation results showed that in the VFL model, compared to the VL model, sarcomere length changed during the isovolumic phase, leading to a decrease in contractile force and shortening of IVRT, which may contribute to improved hemodynamic efficiency. These results indicate that sarcomere length dynamics can mechanically regulate force decay during isovolumic relaxation, even under constant left ventricular volume. This study provides a theoretical framework for understanding the contributions of different layers within the LV wall to diastolic function during the isovolumic relaxation phase.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continuously evolved since its emergence in the human population in 2019. As of 1st August 2025, more than 1,700 Omicron subvariants have been designated by the Pango nomenclature system. The Pango nomenclature system designates a new lineage based on genetic and epidemiological information of SARS-CoV-2 strains. However, there is a possibility that strains that have similar genetic backgrounds and the same phenotype are given different Pango lineage names. In this paper, we propose a new algorithm, called FindPart-w, which can identify groups of viral lineages that share the same relative effective reproduction numbers. We introduced a new lineage replacement model, called the constrained RelRe model, which constrains groups of lineages to have the same relative effective reproduction numbers. The FindPart-w algorithm searches the equality constraints that minimise the Akaike Information Criterion of constrained RelRe models. Using hypothetical observation count data created by simulation, we found that the FindPart-w algorithm can identify groups of lineages having the same relative effective reproduction number in a practical computational time. Applying FindPart-w to actual real-world data of time-stamped lineage counts from the United States, we found that the Pango lineage nomenclature system may have given different lineage names to SARS-CoV-2 strains even if they have the same relative effective reproduction number and similar genetic backgrounds. In conclusion, this study showed that viruses that had the same relative effective reproduction number were identifiable from temporal count data of viral sequences. These findings will contribute to the future development of lineage designation systems that consider both genetic backgrounds and transmissibilities of lineages.

In this research, we create a new fractional-order SEIHRD framework to examine how the Nipah virus moves from one species to another (zoonotic spillover) and how it later spreads throughout a community (via contact with one another) or in a hospital or isolation situation (via entering into a hospital or being placed under quarantine). We used the fractional-derivative formulation of the SEIHRD model to demonstrate memory-based effects related to the progression of an infection and also reflect time-distributed effects associated with surveillance and control measures placed on an infected patient. We first demonstrated that the basic epidemiologic properties of the model were consistent by showing that the solutions of the SEIHRD differential equations will always yield positive and bounded solutions within biologically relevant parameter ranges. We then established the well-posedness of this model by transforming the SEIHRD differential equations into an equivalent integral operator and applying various fixed-point arguments to demonstrate that there will always be unique solution(s) to the SEIHRD differential equations. To evaluate the threshold parameter for the transmission of Nipah virus within a given population we calculated the threshold level through the next generation method to determine the expected number of secondary infections from a new or chronically infected host. One of the main contributions of this work is to include an analysis of the robustness of a given solution to all potential perturbations (i.e., Ulam-Hyers and generalized Ulam-Hyers stability). In addition, we provide analytic results guaranteeing that small perturbations due to approximate modeling, numerical approximation (discretization), or the lack of data fidelity will produce controlled deviations in the solutions. To finish this project, we perform a global sensitivity analysis on uncertain coefficients to evaluate their contribution to the uncertainty of each coefficient and to find out the coefficients that most strongly influence major outcome metrics. This will allow us to develop a priority order for prioritizing spillover control (reduction of human contact and/or isolation), contact reduction, and expenditure of resources towards isolation-related interventions. The resulting framework converts fractional epidemic modeling from a descriptive simulation to a replicable method with robustly defined behavior and equal response prediction.

Human immunodeficiency virus (HIV) persistence remains a major barrier to cure due to the existence of long-lived latent reservoirs that evade immune clearance and persist despite combination antiretroviral therapy (ART). Although ART effectively suppresses viral replication, treatment interruption often leads to rapid viral rebound originating from these latent reservoirs. In this study, we develop a deterministic mathematical model describing the in vivo dynamics of HIV infection incorporating uninfected CD4+ T cells, infected cells, latent reservoirs, deep latent reservoirs, and infectious and non-infectious virions, while explicitly accounting for the therapeutic effects of reverse transcriptase inhibitors (RTIs), protease inhibitors (PIs), and Tat transcription inhibitors. Analytical results establish positivity and boundedness of solutions and derive the effective reproduction number Re using the next-generation matrix approach. Stability analysis shows that the virus-free equilibrium is locally asymptotically stable when Re < 1, while viral persistence occurs when Re > 1. Numerical simulations were performed to investigate therapy interactions, viral rebound following treatment interruption, and the impact of drug efficacy on viral set-points and latent reservoir dynamics. To further explore therapy interactions, three-dimensional viral set-point surfaces and heat maps were generated to examine how combinations of infection inhibition, viral production inhibition, and transcriptional inhibition influence viral dynamics. The simulations reveal that Tat inhibition suppresses viral transcription, thereby reducing the transition of infected cells into productive infection and limiting viral propagation when combined with conventional ART mechanisms. The therapy parameter planes further demonstrate that strong transcriptional inhibition promotes the transition of infected cells into deep latency, supporting the emerging block-and-lock strategy for functional HIV cure. In addition, a three-dimensional eradication boundary surface and therapy cube were constructed to identify regions of parameter space where Re < 1, corresponding to successful viral control. These visualizations show that viral eradication is unlikely when therapies act independently but becomes achievable when multiple therapeutic mechanisms act simultaneously. Overall, the results highlight the critical role of transcriptional inhibition through Tat-targeting therapies in complementing existing ART regimens. By simultaneously suppressing viral replication and promoting deep latency, Tat-based combination strategies may significantly reduce viral rebound and contribute to long-term functional control of HIV infection.

Against the backdrop of global climate change and accelerating population mobility in 2025, chikungunya fever (CHIKF) exhibited a trend of worldwide spread, significantly increasing the difficulty of controlling tropical mosquito-borne diseases. To enhance the precision of intervention strategies, this study developed an age- and sex-structured human-mosquito interaction dynamic model based on data from the largest CHIKF outbreak ever recorded in China, and conducted a targeted analysis of prevention and control strategies. By decomposing the basic reproduction number and examining population heterogeneity, asymptomatic males aged 15-59 years were identified as the core transmission group. Optimal control analysis revealed that the synergistic implementation of three measures-- reducing the effective human-to-mosquito transmission rate, reducing the effective mosquito-to-human transmission rate, and suppressing mosquito population density--could reduce the overall infection rate by 95.7586%. Among these, mosquito population suppression should be prioritized as a universal core strategy; however, its protective effect on females aged 60 years and above was relatively weak, warranting particular attention. The study further demonstrated that asymmetric intensity combinations targeting these three intervention pathways--such as intensity profiles of "10%, 90%, 90%" or "60%, 80%, 90%"--could achieve effective outbreak control. This research elucidates population-specific transmission patterns and key pathways for intervention intensity, providing a theoretical and strategic foundation for the precise control of mosquito-borne diseases. It also provides actionable operational insights to support rapid response and strategy optimization for future emerging outbreaks. Author summaryCHIKF is a mosquito-borne viral disease that is gradually spreading from tropical regions to other areas. To achieve more precise control of this disease, we developed an age- and sex-structured analytical model based on the largest CHIKF outbreak in China, aiming to provide a scientific basis for responding to potential future outbreaks with inherent uncertainties. The study found that asymptomatic males aged 15-59 years were the primary drivers of transmission and should be prioritized as a key population for reducing viral spread in prevention efforts. When evaluating the effectiveness of different intervention strategies, females aged 60 years and above were the least affected by the implemented measures, indicating that this group should strengthen personal protection to lower their infection risk. Among all control measures, mosquito suppression was the most effective, suggesting that vector control strategies should be prioritized in future outbreak responses.

Within many bacterial colonies, persister cells exist as a subpopulation that is tolerant to antibiotics and other stressors, yet not genetically distinct from the rest of the colony. A recent study has proposed epigenetic inheritance as a mechanism that leads to the presence of persister cells. We analyze a nonlocal PDE--ODE model introduced in that study to describe the epigenetic inheritance process and establish its mathematical well-posedness, including existence, uniqueness, and nonnegativity of solutions. We identify a sharp parameter threshold delineating extinction from persistence of the colony: below this threshold the washout equilibrium is globally asymptotically stable, while above it a unique positive equilibrium exists and the population is weakly persistent. Notably, this threshold is independent of the internal community structure.

Variability in gene expression among single cells and growing cell populations can arise from the stochastic nature of protein synthesis, which often occurs in random bursts. This study investigates the variability in the expression of a growth-sustaining protein, whose concentration is regulated by a negative feedback loop due to cell growth-induced dilution. We model the distribution of protein concentration using a Chapman-Kolmogorov equation for single cells and a population balance equation for growing cell populations. For single cells, we derive an explicit solution for the protein concentration distribution in state space and represent it as a Bessel function in Laplace space. For growing populations, we find that the distribution satisfies a Heun differential equation with singular boundary conditions. By addressing the central connection problem for the Heun equation, we quantify the population-level protein distribution and determine the Mathusian parameter, which characterizes population growth. This work provides a comprehensive analytical framework for understanding how stochastic protein synthesis impacts gene expression variability and population dynamics.

Bladder cancer presents significant clinical challenges due to its complex immune microenvironment and highly heterogeneous response to treatments. To create accurate, individualized models of disease progression, we first construct a system of Ordinary Differential Equations (ODEs) that captures tumor-immune interactions. We address the challenge of estimating unknown parameters by performing a rigorous comparative analysis of two heuristic optimization methods: Differential Evolution (DE), a robust global optimization algorithm, and Physics-Informed Neural Networks (PINN), a novel machine learning framework that embeds ODE constraints into its loss function. Our findings provide a critical evaluation of the computational efficiency and accuracy of each method for parameterizing biological ODE systems. This study validates the power of hybrid machine learning approaches in mathematical oncology, yielding a data-driven model of bladder cancer progression with direct potential for optimizing personalized treatment strategies. Author summaryBladder cancer remains a major global health threat, characterized by highly unpredictable responses to treatment and a high likelihood of recurrence. To better predict how a patients disease will progress, researchers use mathematical models that simulate the "war" between cancer cells and the immune system. However, these models are only useful if they can be accurately tuned to a specific patients data--a process called parameter estimation. This task is notoriously difficult because clinical data is often sparse and noisy, making it hard to find the right settings for the model. In this study, we developed a novel computational framework that combines a traditional "survival of the fittest" optimization algorithm (Differential Evolution) with Physics-Informed Neural Networks (PINNs), a specialized architecture designed to embed physical constraints directly into the learning process. By "teaching" the AI the underlying biological laws of cancer growth, our hybrid approach can accurately estimate a patients unique disease parameters even when raw data is limited. We validated this method using a "virtual patient" system derived from real-world clinical trials. Our results show that this hybrid approach provides a more robust and reliable way to personalize cancer models, offering a powerful new tool for doctors to simulate and optimize individual treatment plans before they are even administered.

Mass-conserved reaction-diffusion systems are used as mathematical models for various phenomena such as cell polarity. Numerical simulations of this system present transient dynamics in which multiple stripe patterns converge to spatially monotonic patterns. Previous studies indicated that the transient dynamics are driven by a mass conservation law and by variations in the amount of substance contained in each pattern, which we refer to as "pattern flux". However, it is challenging to mathematically investigate these pattern dynamics. In this study, we introduce a reaction-diffusion compartment model to investigate the pattern dynamics in view of the conservation law and the pattern flux. This model is defined on multiple intervals (compartments), and diffusive couplings are imposed on each boundary of the compartments. Corresponding to the transient dynamics in the original system, we consider the dynamics around stripe patterns in the compartment model. We derive ordinary differential equations describing the pattern dynamics of the compartment model and analyze the existence and stability of equilibria for the reduced ODE with respect to the boundary parameters. For a specific parameter setting, we obtained results consistent with previous studies. Moreover, we present that the stripe patterns in the compartment model are potentially stabilized by changing the parameter, which is not observed in the original system. We expect that the methodology developed in this paper is extendable to various directions, such as membrane-induced pattern control.

The placenta is a fundamental organ for human reproduction, facilitating fetal growth via the exchange of oxygen, nutrients and waste products between mother and fetus, by means of a dense network of fetal villi bathed in maternal blood that flows through the intervillous space (IVS). However, despite its role in adverse pregnancy outcomes associated with impaired maternal-fetal transfer, the influence of placental structure on maternal haemodynamics and the delivery of oxygen and nutrients to the developing fetus is not well understood. This study employs computational fluid dynamics within physiologically informed 2D and 3D whole-organ-scale placental geometries, whose features are informed by recent ex vivo experimental and {micro}CT data, to examine comprehensively the influence of placental anatomy on maternal flow and transport in the IVS. In particular, we consider in detail the impact of the number and placement of maternal decidual arteries and veins that supply and drain the IVS, the location and height of so-called septal walls that loosely separate the placenta into functional units (cotyledons) and sub-units (lobules), and the density of the fetal villous trees as reflected in the rate of uptake of dissolved solutes from the maternal blood and the resistance to flow. We first exploit the computational efficiency of simulation in a representative 2D geometry to study in detail the sensitivity of haemodynamic markers to these parameters. These results guide our 3D study which reveals that the flow, transport and oxygen uptake are strongly influenced by placental structure, and exposes the vein-to-artery ratio as a key indicator of placental efficiency, regulating a trade-off between a preferential maternal flow environment and fetal oxygen uptake. Conversely, the location and height of the septal walls, a feature that is not well studied, have minimal systematic impact on the macroscopic haemodynamic and transport measures considered here. We also introduce a reduced model, for which analytical progress can be made, and demonstrate its utility in exposing key drivers of maternal-fetal transport. Author summaryIn our study, we explored how placental structure affects maternal blood flow and the delivery of oxygen to the fetus. The placenta is a complex and vital organ, but we do not fully understand how its morphology influences its functional efficiency. This is a critical gap in our knowledge as placental blood flow is implicated in serious pregnancy complications such as fetal growth restriction and pre-eclampsia. To investigate the effect of structure, we used detailed 2D and 3D models of the placenta. These models allowed us to simulate maternal blood flow and oxygen transport therein while changing various structural features, such as the number and location of placental veins and the density of the intervillous space. We found that the ratio of veins to arteries is a key factor in determining the maternal blood flow speed and the amount of oxygen the fetus receives. A higher number of veins relative to arteries helps create the slow-flow environment necessary for a healthy pregnancy, but an excessive ratio can reduce the efficiency of oxygen uptake. Our findings suggest that the vein-to-artery ratio could be an important indicator of placental health, offering a new perspective for understanding and enabling effective early intervention treatments of pregnancy complications.

Visceral leishmaniasis (VL) is considerably more severe among individuals infected with human immunodeficiency virus (HIV), leading to higher parasite loads, frequent relapse, and increased mortality. To examine the epidemiological interaction between the two diseases, we develop a comprehensive VL-HIV co-infection model that incorporates transmission pathways, treatment effects, and relapse dynamics. The model is parameterized using real-time data from Bihar, India, including monthly VL-only and VL-HIV co-infected cases and annual HIV prevalence data. Our analysis shows that HIV infection drives the resurgence and persistence of VL even in settings where VL alone would not sustain transmission, underscoring the amplifying effect of HIV-induced immunosuppression on VL dynamics. We further demonstrate that increasing HIV treatment coverage substantially reduces co-infection prevalence and lowers VL relapse rates. Numerical simulations and optimal control analysis highlight the effectiveness of integrated intervention strategies that combine awareness, treatment enhancement, and vector control. Overall, this study emphasizes the need for coordinated VL and HIV control programs and provides data-driven guidance for designing sustainable intervention strategies in endemic regions.

With significant population fractions in many societies who refuse vaccines, it is important to reconsider how vaccination is incorporated into compartmental epidemiology models. It is still most common to apply the vaccination rate to the entire class of susceptibles, rather than to use the more realistic assumption that the vaccination rate function should depend only on the population of susceptibles who are willing and able to receive a vaccination. This study uses a simple generic disease model to address two questions: (1) How much error is introduced in key model outcomes by neglecting vaccine unwillingness?, and (2) Can the error be reduced by incorporating vaccine unwillingness into the vaccination rate constant rather than the rate diagram? The answers depend greatly on the time scale of interest. For the endemic time scale, where longterm behavior is studied with equilibrium point analysis, the error in neglecting unwillingess is large and cannot be improved upon by decreasing the vaccination rate constant. For the epidemic time scale, where the first big epidemic wave is studied with numerical simulations, the error can still be significant, particularly for diseases that are relatively less infectious and vaccination programs that are relatively slow.

Infectious disease surveillance systems in tropical countries show that respiratory disease incidence generally manifests as year-round activity with weak fluctuations and irregular seasonality. Previously, using a ten-year time series of influenza-like illness (ILI) collected from outpatient clinics in Ho Chi Minh City (HCMC), Vietnam, we found a combination of nonannual and annual signals driving these dynamics, but with unknown mechanisms. In this study, we use seven stochastic dynamical models incorporating humidity, temperature, and school term to investigate plausible mechanisms behind these annual and nonannual incidence trends. We use iterated filtering to fit the models and evaluate the models by comparing how well they replicate the combination of annual and nonannual signals. We find that a model including specific humidity, temperature, and school term best fits our observed data from HCMC and partially reproduces the irregular seasonality. The estimated effects from specific humidity and temperature on transmission are nonlinearly negative but weak. School dismissal is associated with decreased transmission, but also with low magnitude. Under these weak external drivers, we hypothesize that stochasticity makes a strong sub-annual cycle more likely to be observed in ILI disease dynamics. Our study shows a possible mechanism for respiratory disease dynamics in the tropics. When the external drivers are weak, the seasonality of respiratory disease dynamics is prone to the influence of stochasticity.

Common knowledge states that the spontaneous somatic evolution of a normal tissue may lead to a tumor. Once the tumor is formed, it naturally evolves towards a state of higher malignancy. On the other hand, perfect gene expression markers for normal tissue and tumor--the so-called N-genes and T-genes--were recently introduced. We join these two pieces of knowledge in order to argue that: 1) Only N-markers participate in the spontaneous dynamics of a normal tissue. The number of active markers decreases as the tissue approaches the transition point where it becomes a tumor. 2) Only T-markers participate in the spontaneous dynamics of tumors. The number of markers increases as the tumor becomes more malignant. 3) Both sets of genes are connected by the so-called NT-genes, i.e., genes that are simultaneously N- and T-markers. They should play a crucial role at the transition point and, possibly, when the tumor is exposed to a drug or therapy. 4) The pathways or mechanisms protecting the normal tissue from becoming a tumor may be described by a small perfect panel of N-genes. 5) The pathways or mechanisms guiding the evolution of tumors in a tissue may be described by a small perfect panel of T-genes. We illustrate the above statements with the analysis of expression data for prostate adenocarcinoma, one of the most heterogeneous tumors. In this case, there are about 1000 N-genes and 6000 T-genes, and the perfect N- and T-panels contain 11 and 8 genes, respectively. Additionally, we provide examples from lung adenocarcinoma and liver hepatocarcinoma.

Monkeypox viral disease has been and continues to be a global public health concern. Currently, there are existing, though minimal measures to manage mpox and any future outbreaks. Relying on data-driven modeling for early detection of mpox and prediction of possible cases and deaths in the presence of an outbreak is thus imperative. The present study forecasted global mpox virus cases and deaths in Asia, Africa, Australia, Europe, North America, Oceania, and South America. Three forecasting models (deep neural network, gradient boosting, and polynomial regression) were trained on data from the seven geographical regions. The performance of the three models was assessed using coefficient of determination, mean squared error, root mean squared error, and mean absolute scaled error across each region. Prediction using the deep neural network revealed a potential of higher mpox deaths in Africa and higher mpox cases in South America. Prediction using gradient boosting showed a potential of mpox deaths in Africa and higher mpox cases in Asia and North America. Prediction using polynomial regression revealed a potential of higher mpox deaths in Africa and Asia while rapid rises in mpox cases from 2025 to 2028 were anticipated in all regions except Asia in case of a monkeypox outbreak. For the three models, the tree-based ML model (gradient boosting) outperformed the statistical model and deep learning model by R2 and MSE in predicting mpox case counts across all the seven geographical regions. This study showcases the worth in using data-driven modelling to predict emerging and re-emerging infectious diseases such as mpox.

Cardiovascular disease is one of the leading causes of death worldwide, with myocardial infarction (MI) being a major cause of both morbidity and mortality among cardiovascular patients. MI Patients face a higher risk of cardiovascular disease recurrence afterwards. Therefore, accurately predicting the risk of recurrence and identifying key risk factors are crucial for clinical decision-making. In this paper, we consider the interrelationships among cardiovascular factors from a systemic perspective. We first construct a differential network for each patient to capture individual-specific deviations in factor relationships and propose a novel method, termed Causal Factor-aware Graph Neural Network (CFGNN), which integrates factor interactions to predict the recurrence risk of MI patients while uncovering key risk factors from a causal perspective. Experimental results demonstrate that CFGNN performs well on hospital-derived datasets in real world, effectively identifying several key risk factors. This method not only deepens our understanding of cardiovascular disease, but also paves the way for more targeted and effective interventions.

Weeds present a major challenge to agricultural productivity, competing with crops for critical resources like water, nutrients, and sunlight, resulting in significant yield reductions. Prompt weed identification is essential for enabling effective control strategies, such as the application of herbicides or mechanical removal, to minimize their impact on crop growth. This research focuses on developing a deep learning approach based on game theory for detecting weeds. Using CWFID dataset captured at various times and days, along with multispectral data in the visible and near-infrared spectrum, the study aims to improve early detection methods for more efficient weed management in agricultural settings. A novel segmentation technique for weed regions is introduced, employing a zero-sum game theory model to reconcile conflicting classifications from different weed detectors. These regions are treated as zones of conflict between weeds and crops, with each detector representing a different strategy. By defining an appropriate utility function, the method identifies the Nash equilibrium, effectively minimizing false positive detections of weeds.

Mainland China experienced multiple waves of COVID19 pandemic during 2020 2022, driven by emerging variants and changes in public health and social measures (PHSMs). We developed a hypergraph-based Susceptible Vaccinated Exposed Infectious Recovered Susceptible (SVEIRS) model to reconstruct epidemic dynamics across 31 provinces, capturing transmission heterogeneity associated with clustered contacts. We assessed key characteristics of transmission at national and provincial levels during four outbreak periods: initial, localized predelta, Delta, and widespread Omicron, which accounted for 96.7% of all infections. We found significant diversity in transmission contributions across cluster sizes, with a small fraction of larger clusters responsible for a disproportionate share of infections. Counterfactual analyses showed that reducing clustersize heterogeneity, while holding overall exposure constant, could have lowered national infections by 11.70 to 30.79%, with the largest effects during Omicron period. Ascertainment rates increased over time but remained spatially heterogeneous with a range: (14.40, 71.93)%. Population susceptibility declined following mass vaccination (to 42.49% in Aug 2021, nationally) and rebounded (to 89.89% in Nov 2022) due to waning immunity with variations across the provinces. Effective reproduction numbers displayed marked temporal and spatial variability, with higher estimates during Omicron. Overall, these results highlight critical role of group contact heterogeneity in shaping epidemic dynamics.